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The cosmic first order phase transition from quarks to hadrons, occurring a few microseconds after the Big Bang, 
would lead to the formation of quark nuggets which would be stable on a cosmological time scale, if the associated 
VjQ ' baryon number is larger than a critical value. We examine the possibility that these surviving quark nuggets may not 
04 . only be viable candidates for cold dark matter but even close the universe. 
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t • . According to the wisdom of the standard model, the universe underwent, a few microseconds after the big bang, a 

phase transition from quarks to hadrons. The cosmological implications of this, presumably first order, phase transition 
would be far-reaching. Schramm and collaborators (l],|j (see also ||) argued that the associated fluctuations may lead 
fS to the formation of primordial black holes, which could be as large as Mq, the solar mass, for fluctuations around 
the horizon scale. They recently suggested Q that these black holes could even be the candidates for the Massive 
Compact Halo Objects (MACHO 's) which had of late been discovered in the halo of the Milky way, in the 

direction of the Large Magellanic Cloud (LMC) through their gravitational microlensing properties. On the other 
C$ • hand, a first order phase transition scenario involving bubble nucleation at a critical temperature T c ~ 100 — 200 MeV 
should lead to the formation of quark nuggets (QN) Q, made of u, d and s quarks at a density > nuclear density. If 
these primordial QN's existed till the present epoch, they could be possible candidates for the dark matter |J. Such 
a possibility would be aesthetically rather pleasing, as it would not require any exotic physics nor would the success 
of the primordial (Big Bang) nucleosynthesis scenario be affected |7|-|l0|. 

The central question in this context then is whether the primordial QN's can be stable on a cosmological time 
scale. In a recent work |fl| using the chromoelectric flux tube model, we have demonstrated that the QN's will 
survive against baryon evaporation, if the baryon number of the quark matter inside the nugget is larger than 10 42 . 
For reasons explained in [ pr| , this estimate is rather conservative. Sumiyoshi and Kajino ]l3] have estimated that 
a QN with an initial baryon number ~ 10 39 would survive against baryon evaporation. It should be noted at this 
point that the horizon limit on the baryon number at that primordial epoch is around 10 49 . The conclusion from 
these calculations is that the larger primordial QN's within the baryon number window 10 39-40 < Nb < 10 49 are 
indeed cosmologically stable. (Also noteworthy in this context is the observation that these numbers are sufficiently 
above the limit (Nb ~ 10 21 ) which was obtained by Madsen some time ago, below which the QN's could interfere 
with the outcome of Big Bang nucleosynthesis.) It is therefore most relevant to ask what fraction of the dark matter 
could be accounted for by the QN's. The central issue that we wish to address in this letter is whether we can have 
a scenario where the universe would be closed with QN's of baryon number within the above range. In other words, 
can the proverbial cosmological dark matter, containing 90% or more of all the matter in the universe, be made up 
entirely of such QN's ? 

It is well known that in a first order phase transition, the quark and the hadron phases co-exist. This configuration 
would be referred to in the following as the mixed phase. In the quark phase, the universe consists of leptons, photons 
as well as the quantum chromodynamic (QCD) degrees of freedom (massless quarks, anti-quarks and gluons) and 
is described by, say, the MIT bag equation of state with an effective degeneracy <?q w 51.25. Note that the baryon 
number in this phase is carried entirely by the quarks. The hadronic phase contains baryons, mesons, photons and 
leptons and is described by an equation of state corresponding to massless particles with an effective degeneracy 
9h = 17.25. 

The evolution of the universe in the mixed phase at the critical temperature T c of the phase transition is governed 
by the Einstein equation in the Robertson- Walker spacetime, as described below : 
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where e is the energy density, P the pressure and G the gravitational constant. Combining eqs. (1) and (2) with the 
equation of state, we determine the scale factor R(t) and the volume fraction of the quark matter f(t) in the mixed 
phase as 
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3m^/87r_B is the characteristic time scale for the QCD phase transition in the early universe 



where r = g q /g h , t c - y u„* pJ 

and t{ is the time when phase transition starts. The bag constant, B 



(245) 4 MeV 4 in our calculations. 
In the mixed phase, the temperature of the universe remains constant at T c , the cooling due to expansion being 
compensated by the liberation of the latent heat. In the usual picture of bubble nucleation in first order phase 
transitions, hadronic matter starts appearing in the quark matter as individual bubbles. With the progress of time, 
more and more hadronic bubbles form, coalesce and eventually percolate to form a network of hadronic matter which 
traps the quark phase into finite domains (l3| . The time when the percolation takes place is called the percolation 
time t p , determined by a critical volume fraction f c , (f c = f(t p )) of the quark phase. 

Detailed numerical studies on percolating systems yield the result that for bubbles with the same radial size, f c 
is ~ 0.3 |L4|,|l5|. We would also use the same value of f c here. For the sake of simplicity, we would also assume that 
the trapped quark domains are all of the same size. (In general, this is not true; there ought to be a distribution 
of domain sizes. This is an involved problem which we plan to look into in future. The conclusions of this work, 
however, are not expected to be drastically altered.) From eq.(4), we get t p — tj = 0.08i c . The values of various time 
scales in the present case are t c = 46 /isec and t p = 27 fxsec. 

In an ideal first order phase transition, the fraction of the high temperature phase decreases from the critical 
value f c , as these domains shrink. For the QCD phase transition, however, these domains should become QN's 0] 
and as such, we may assume that the lifetime of the mixed phase tf (i.e., the time when the cooling due to expansion 
starts to dominate again and the temperature of the universe starts falling), is ~ t p . 

The probability of finding a domain of trapped quark matter of co-ordinate radius X at time t p is given by p3[ , 



P(z, x p ) = exp 
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where z = XR(ti)/vt c , x — t/t c , r(x) — R{x) / R{xi) and I(x) is the rate of nucleation per unit volume, v is the radial 
growth velocity of the nucleating bubbles, which we left as a parameter. y(x p , x) is given by the following equation 
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Various authors 1 13 16|-|l8[] have proposed different nucleation rates for the cosmic QCD phase transition. Let us 
start with the prescription of ref . fl3|| , where the nucleation rate is given by the following expression 

J(t) = r T S(t - U) (7) 

where the prefactor of the thermal nucleation rate is determined from the requirement P(0, t p ) = f c = 0.3, yielding 
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with 



1.2. This leads to 



P(z, x p ) = exp 



(9) 
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Fig.l: Probability of finding a domain of co-ordinate size X in the quark phase at time t as a function of RiX/vt c 
with the nucleation rate given by eq.(7). 

In fig.l, we plot the probability P(z, x p ) as a function of z. The radius of the trapped quark domain is determined 
by the length scale where the probability falls to f c /e. This implies 
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The number density of the QN's (uqn) can be obtained by using the relation uqnVqn — fc as 
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In an idealised situation where the universe is closed by the baryonic dark matter trapped in the QN's, we should 
have, 



N§(t p )=N2 N n QN V H (t p ) 



(12) 



where iVjf (t p ) is the total number of baryon required to close the universe (fig = 1) at t p , N® N is the total number 
of baryons contained in a single quark nugget and Vjj(i p ) = 47r(ci p ) 3 /3 is the horizon volume. 
Demanding that v/c < 1/ \/3, we get 
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Since the usual baryons constitute only ~ 10% of the closure density, a total baryon number of 10 within the 
horizon at a temperature of ~ 100 MeV would close the universe baryonically. This would require N to be < 10 45 ' 3 , 
which is within the survivability limit of QN's mentioned earlier. 

To study the sensitivity of our results to the nucleation rate, we evaluate the density of QN's for the different 
nucleation scenarios referred to above. 

In the nucleation scenario of Cottingham et al |L6j the density of QN's is found to be larger than the previous 
case n QN = 42947/(wt p ) 3 and the corresponding limit on N§ N is given by N$ N = lO^TVf (t p ). For N§ (t p ) - 10 50 , 
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to be < 10 . Although we have used the exact nucleation rate obtained in eq.(8) of |]16j, it is straightforward 



to show that it can also be approximated to a delta function. With the nucleation scenario of refs. |L7| and 
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< 10 and 10 respectively for the value of the surface tension = 50 MeV/fm 
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It thus appears that the upper limit on the baryon number of QN's that would close the universe baryonically 
is not very sensitive to the nucleation mechanism ( little about which is confidently known due to the rudimentary 
knowledge of the dynamics of the QCD phase transition) and all available estimates point to the real possibility that 
stable QN's could not only be a possible candidate for cold dark matter but they could even close the universe with 
total baryon numbers within the window mentioned earlier. 



For the sake of completeness, we would like to mention here that QN's with baryon number lower than the 
survivability window would evaporate rather quickly, leaving a large baryon inhomogeneity. This would still not 
pose much problem for the cosmological scenario, as this overdensity would dissipate, primarily due to neutrino 
inflation upto temperatures ~ 1 MeV and then by baryon diffusion, which becomes the dominant mechanism at lower 
temperatures Jf9|-pl[. Even if the initial overdensity due to the evaporating QN's could be as high as 10 10 ~ 12 , it is 
found to go down by several orders of magnitude by the time nucleosynthesis starts. These details would be reported 
elsewhere plj . 

In order for the stable QN's to be a viable candidate for the cosmological dark matter, they must gravitationally 
clump, just like normal matter. It may however be asked what, if any, observational signatures could these stable 
QN's have at the present time. To this end, we may estimate the rate at which such QN's may collide with the earth, 
in the extreme case where all these QN's were floating around in the universe. For the various scenarios [^3 
considered here, the rate of such collisions (density of QN'sx velocity of QN's (vqn ~ 2 x 10 7 cm/sec) xsurface area 
of the earth (~ 8 x 10 17 cm 2 )) turns out to be one in ~ 10 15 — 10 16 years for v = c/v3. The collision rate could be 
much larger for smaller velocity of the bubble growth as the rate varies as 1/v 3 . 

More realistically, however, the scenarios proposed here would get substantial experimental support, albeit some- 
what indirect, if a sizable flux of stable strangelets ( heavy nuclear objects with very abnormal charge-to-mass ratios 
) could be detected in the cosmic rays. We may recall that there has been a report in the literature |23] that such 
objects may indeed form part of the extragalactic cosmic rays. Extensive plans are under way at this time to look 
for such objects in high altitude experiments |2|,Q, the outcome of which would go a long way to shed light on the 
nature of dark matter. 

We are grateful to Dr. Pijushpani Bhattacharjee for stimulating discussions. 
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